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Second and Third Virial Coefficients for Hydrogen 
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Second and third virial coefficients for parahydrogcn have been derived from closely 
spaced PVT data from 24 to 100 °K. They are in good agreement at 100 °K with published 
data for normal hydrogen. Analytical representations of the combined data from about 
20 to 423 °K are presented which may be useful in computation of thermodynamic functions 
of the gas. These formulas are related to those resulting from the use of the Lennard-Jones 
potential. 



1. Introduction 

As part of an extensive program for determining 
the thermodynamic and transport properties of fluid 
parahydrogen, we have recently published new data 
on the P-V-T surface from 15 to 100 °K [I]. 1 Ear- 
lier, we gave preliminary values of the second and 
third virial coefficients, B(T) and C(T), in a paper 
describing the apparatus [2]. 

In this paper we present final values of the virial 
coefficients derived from the smoothed compressibil- 
ity data [1]. These differ but slightly from the pre- 
liminary values. Also we present and examine ana- 
lytical representations of B(T) and (\T) which may 
be useful for computing thermodynamic functions 
of the gas. 

In a forthcoming paper [3] we have used the fol- 
lowing arbitrary representations [4] to fit our closely 
spaced data in order to facilitate computations of 
thermofunctions below 100 °K: 



B=JZ b t T 



f^JZdT 



On the other hand, in this paper use is made of other 
published virial coefficients as well as those from this 
laboratory to provide a more extensive tabulation 
from about 20 to 423 °K to which are fitted expres- 
sions having forms suggested by the use of the 
Lennard-Jones potential. All coefficients were de- 
termined by the method of least squares. 



2. Derivation of the Virial Coefficients 

The virial expansion [5] may be rearranged to a 
form which, when truncated, is linear in density on 
isotherms, convenient for graphical or analytical de- 



1 Figures in brackets indicate the literature references at the end of this paper. 



termination of the second and third coefficients, 



$=(Pv-rt)v=rtb+rtc/v+ 



(i) 



Parameters RTB and RTC of eq (1) as truncated 
have been determined on isotherms by least-squares 
using the reported compressibility data [1] and the 
orthobaric densities [6]. The number of datum 
pairs, n, on each isotherm is given in column (a) of 
table 1. It was found by trial that data for densities 
above 0.007 g mol/cm 3 diverge from the linear be- 
havior of eq (1). These data were omitted, with the 
result that a maximum of eight points per isotherm 
were applicable above 29 °K. (The above maximum 
density for validity of eq (1) with only two terms is 
less than half the critical density of 0.0156 g mol/cm 3 
[6].) Derived values of B and fare given in tables 
2 and \\, respectively, at integral temperatures from 
24 through 100 °K. 

Equation (1) also was used for smoothing and 
interpolation of compressibility data to the highest 
densities by admitting as many terms as required by 
the experimental precision [1]. At temperatures 
below critical (32.9S °K), a single isotherm of eq (1) 
was used to represent both vapor and compressed 
liquid. The number of datum pairs, n' ', and number 
of terms, N, in polynomial (1) are given in column (b) 
of table 1 which includes also the maximum density 
multiplied by 10 3 for each isotherm. With increas- 
ing number of terms, the mean deviation decreases 
more at an odd number than at an even number of 
terms. Since coefficients of the higher-order terms 
of these divergent, alternating polynomials show no 
regular temperature-dependence, the expansions can- 
not be interpreted as virial equations [7]. Neverthe- 
less, values of B and C obtained in this way are pre- 
sented in appendix, table 1 , to meet the commonplace 
question of the effect upon them of an increased num- 
ber of terms in eq (1). it may be seen that values 
of B are essentially the same as in table 2. Values 
of (\ while of similar magnitude to those in table 3, 
scatter badly. 
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Table 1. Number of datum pairs, n, and of terms N, for eq(l) 





( a ) 




00 






(a) 




(b) 




T,°K 










T, °K 












n 


n' 


N 


103Mnin 


n 


n' 


N 


10 3 /tfmin 


24 


3 


24 


9 


44.1 


40 


8 


45 


11 


40.4 












42 


8 


42 


11 


40.3 


25 


3 


25 


9 


44.1 


44 


8 


41 


9 


39.7 


26 


4 


25 


9 


43.8 


46 


8 


40 


9 


39.1 


27 


5 


27 


9 


43.8 


48 


8 


39 


9 


38.4 


28 


6 


29 


9 


43.4 












29 


7 


31 


9 


43.4 


50 


8 


38 


9 


37.6 












55 


8 


37 


9 


36.6 


30 


8 


33 


11 


42.9 


60 


8 


36 


9 


35.9 


31 


8 


35 


13 


42.9 


65 


8 


34 


9 


34.6 


32 


8 


38 


13 


42.5 


70 


8 


33 


7 


34.0 


33 


8 


47 


15 


42.4 












34 


8 


48 


13 


42.0 


75 


8 


30 


7 


31.8 












80 


8 


29 


7 


30.8 


35 


8 


48 


11 


42.0 


85 


8 


28 


7 


29.8 


36 


8 


47 


11 


41.5 


90 


8 


27 


7 


28.9 


37 


8 


47 


11 


41.5 


95 


8 


26 


5 


28.1 


38 


8 


46 


11 


40.9 












39 


8 


46 


11 


40.9 


100 


8 


24 


5 


26.6 



* For N=2 and 0.001 <l//;<0.007 g mol/cm*. 
b To the highest experimental densities, KP/r* 



Table 2. Derived and calculated values of 


B, cm z lg mol 






Calculated 






Calculated 


T,°K 


Derived 






T, °K 


Derived 








(2a) 


(2b) 


(2a) 


(2b) 


15 




-219.4 


-222. 3 


50 


-33. 39 


-33. 25 


-33. 47 


16 


-204. 2 


-200.8 


-202. 8 


55 


-27.48 


-27.29 


-27. 53 


17 




-184.7 


-186.1 


60 


-22.70 


-22.43 


-22.67 


18 


-172.9 


-170.6 


-171.5 


65 


-18.64 


-18.40 


-18.63 


19 




-158.2 


-158. 7 


70 


-15.22 


-15.01 


-15.23 


20 


-148. 8 


-147.1 


- 147. 4 


75 


-12.42 


-12.13 


-12.32 


21 




-137.2 


-137. 4 


80 


-9.88 


-9.65 


-9.82 


22 


-129.7 


-128.4 


-128.4 


85 


-7.63 


-7.50 


-7.63 


23 




-120.3 


-120.3 


90 


-5.66 


-5.61 


-5.72 


24 


-112.8 


-113.1 


-113.0 


95 


-3.99 


-3. 95 


-4.02 


25 


-106. 2 


-106. 5 


-106.3 


98.15 


-3.06 


-2.99 


-3. 05 


26 


-100.3 


-100.4 


-100. 3 


100 


-2. 52 


-2.47 


-2.51 


27 


-94. 80 


-94. 88 


-94. 74 


103. 15 


-1.69 


-1.62 


-1.65 


28 


-89. 66 


-89. 78 


-89. 65 


113. 15 


+0.67 


+0.73 


+0.75 


29 


-85. 03 


-85.08 


-84. 96 


123. 15 


2.63 


2.65 


2.73 


30 


-80. 73 


-80. 72 


-80. 62 


138. 15 


5.01 


4.96 


5.09 


31 


-76. 75 


-76.68 


-76.60 


153. 15 


6.89 


6.76 


6.94 


32 


-72. 99 


-72.93 


-72.86 


173. 15 


8.84 


8.62 


8.85 


33 


-69. 53 


-69.43 


-69. 39 


198. 15 


10.65 


10.37 


10.62 


34 


-66. 22 


-66. 16 


-66. 14 


223. 15 


11.98 


11.68 


11.93 


35 


-63. 17 


-63. 09 


-63. 10 


248. 15 


12.97 


12.70 


12.94 


36 


-60. 26 


-60. 22 


-60. 25 


273. 15 


13.76 


13.51 


13.72 


37 


-57. 54 


-57. 53 


-57. 58 


298. 15 


14.38 


14.17 


14.34 


38 


-54. 99 


-54. 99 


-55.06 


323. 15 


14.87 


14.71 


14.85 


39 


-52.60 


-52.60 


-52. 69 


348. 15 


15.27 


15.17 


15.25 


40 


-50. 32 


-50. 34 


-50. 45 


373. 15 


15. 60 


15.56 


15.59 


42 


-46. 19 


-46. 19 


-46. 33 


398. 15 


15.86 


15.90 


15.87 


44 


-42. 50 


-42. 46 


-42. 62 


423. 15 


16.08 


16.19 


16.10 


46 


-39. 18 


-39.09 


-39. 28 










48 


-36. 17 


-36. 03 


-36. 24 











Uncertainty in B ) estimated from eq (1) for low 
densities, is 

dB=[2Z-l + C/v 2 ] 8v-vZ8T/T (la) 

where Z=Pv/RT, or approximately 

8B~v[8v/v-5T/T]. (lb) 

Experimental uncertainty in v was estimated to be 
0.2 percent at low densities [2]. Examination of the 
sensitivity of the shape of isotherms of (Pv-RT)v 
to errors in density, however, indicates a precision 
near 0.02 percent. This is illustrated by figure 1 



Table 3. Derived and calculated third virial coefficients 


T°K 


C. (cmVg mol) 2 


T,°K 


C.Ccmtyg mo 1)2 




Derived 


Calc'd. 




Derived 


Calc'd. 


20 




-405 

+218 
680 
1018 
1259 


50 
55 
60 
65 
70 


964 
889 
838 
785 
743 


968 


21 




893 


92 




835 


23 




789 


24 


1207 


750 


25 
26 

27 
28 
29 


1402 
1580 
1627 
1612 
1615 


1425 
1533 
1596 
1625 
1629 


75 
80 
85 
90 
95 


726 
694 
659 
636 
624 


718 
690 
665 
643 
624 


30 
31 
32 
33 
34 


1600 
1585 
1550 
1516 
1466 


1614 
1586 
1549 
1507 
1463 


98.15 
100 

103. 15 
113.15 
123. 15 


530 
609 
560 
540 
560 


613 
606 
596 
567 
542 


35 
36 
37 
38 
39 


1426 
1377 
1331 
1290 
1252 


1418 
1373 
1330 
1290 
1251 


138. 15 
153. 15 
173. 15 
198. 15 
223. 15 


540 
522 
500 
458 
437 


510 
483 
454 
424 
399 


40 
42 
44 

46 
48 


1209 
1144 
1091 
1046 
1005 


1215 
1151 
1095 
1047 
1005 


248. 15 
273. 15 
298. 15 
323. 15 
348. 15 


415 
404 
370 
340 
313 


378 
360 
345 
331 
319 








373. 15 
398. 15 
423.15 


303 
310 
302 


308 
298 
289 



in which the upper curve is an experimental linear 
isotherm, and the lower curve shows the effect of an 
artificially introduced error of 0.263 percent in 
density [2]. As a further illustration figure 2 gives 
deviations of the data from the least-squared, 
truncated virial expression at 44 °K and compares 
these with the boundary for a deviation of ±0.02 
percent in density. From such considerations, the 
precision of these lowest experimental densities 
appears to be within 0.03 percent. Since this applies 
to all isotherms, the corresponding error in B should 
be roughly independent of temperature. For the 
lowest density, #=1000 cm 3 /g mol; uncertainty 
in B thus is estimated to be about 0.3 cm 3 /g mol. 
Systematic deviations in the temperature-depend- 
ence of B, on the other hand, may be sought in 
absolute deviations of the NBS temperature scale 
for platinum resistance. For ST=0.02 °K, ?;=1000 
cm 3 /g mol, and 25° < 2X100 °K, the range of 
uncertainty is 0.8<55<0.2 cm 3 /g mol. 

Since B is derived as the intercept of linear plots of 
(Pv/RT—l)v versus 1/v, it may be seen that all of 
the above arguments for precision apply equally 
to absolute errors in v or in T. Absolute uncertainty 
in B therefore is estimated to be in the range 1.1 
cm 3 /g mol at 25 °K to 0.5 cm 3 /g mol at 100 °K. 

To the above virial coefficients for parahydrogen 
from 24 to 100° in tables 2 and 3 have been added 
values for normal hydrogen at higher and lower 
temperatures in order to provide a more extensive 
set. Of the published values of virial coefficients for 
normal hydrogen [4, 8, 9, 10, 11, 12], values from 
Woolley, Scott, and Brickwedde [9] below 24 °K 
are included in table 2, while values from Michels, 
de Graaff, and Ten Seldam [10] at nonintegral 
temperatures above 98 °K are given in tables 2 and 
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Figure 1. Effect of systematic error in density [2] upon the 
linear behavior of eq. (1), for the 83 °K isotherm. 

The notation is Z=Pv/RT, Circles: experimental linear isotherm. Triangles: 
same data with artifically introduced error of 0.263 percent. 
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Figure 2. Deviations of data, <£= (Z-l)RTv, where Z = 
Pv/RT, from eq. (1) at U °K. 
Dashed lines correspond to ±0.02 percent systematic error in density. 



3. Uncertainty of the data of Michels et al., was 
estimated by them to be 0.15 cm 3 /g mol in B ) and 
about 15 percent in C. Excellent agreement of the 
independent data at 100 °K suggests that the virial 
coefficients of these hydrogen modifications may be 
indistinguishable at this and higher temperatures. 
Small differences have been detected at lower 
temperatures [13]. 

3. Representation of Second Virial 
Coefficient 

Figure 3 presents the derived data of table 2 as a 
function of T~ 5/i . The data for normal hydrogen 
at low temperatures are seen to deviate from the 
straight line extrapolated from the parahydrogen 
data, in agreement with the experimental comparison 
of the two modifications made by Beenakker et al. 

[13]. 

Systematic deviations persist with all analytical 
representations which have been investigated for 
the temperature-dependence of B. It therefore does 
not appear possible to select a form which is best for 
all purposes from among the following two, which 
utilize two and four constants, respectively, with 
notation x=T /T: 

B=B [l-x 5/ % T o =109.83 °K, 

# =19.866 cm 3 /g mol, A=0.125 cm 3 /g mol. (2a) 



B=J2 BiX m - 1)n , 7 7 o=109.781 °K, 

^ = +42.464, £ 2 = -37.1172, 

B z = -2.2982, B 4 = -3.0484, 
A=0.066 cm 3 /g mol. 



(2b) 



Mean deviations, given for each equation above, are 
defined as 

^(n-Ny'itlB-B^l 

wherein ^=49 is the number of datum pairs, and N 
is the number of constants in the equation. Calcu- 
lated results from these equations are given in table 
2, and the individual deviations are plotted in 
figure 4. 

The form of (2a) is that used by Keesom for 
helium [14], the constant high-temperature limit 
corresponding to the rigid-sphere model for molecu- 
lar interaction [5]. The value of B in (2a) is com- 
parable with the volume 22.65 cm 3 /g mol of solid 
normal hydrogen at 4.2 °K [9]. The rather good 
fit provided by this simple equation is indicated by 
figure 3. Equation (2b) consists of the leading 
terms of an expansion derived from the Lennard- 
Jones potential [5]. Appendix, table II, gives 
coefficients for that expansion with increasing num- 
ber of terms. (The form of the expansion given in 
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Figure 3. Second virial coefficient, B, versus T~ 5fi . 

Open circles are for paruhydrogen. Filled circles for normal hvdrogen below 21 °K are from [91 and above 100 °K from 
[10]. 
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Figure 4. Deviations of second virial coefficient, B, in cm d /g mol, from eq (2a), open circles 

and from (2b), filled circles. 



the appendix does not include the root, T .) It is 
seen that four terms give an acceptable representa- 
tion. Whereas use of additional terms improves 
the fit, the similarity in form of the two sets of devia- 
tions in figure 4 suggests the presence of systematic 
experimental errors. If such errors are present, it 
is doubtful that use of higher terms would be justi- 



fied. The form of the deviations in figure 4 having 
been found also with analytical representations of 
B other than polynomials, it is improbable that they 
arise from the selected polynomial forms of (2a) and 
(2b). It appears equally improbable that they 
arise entirely from deviations of the NBS tempera- 
ture scale used. 
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4. Representation of Third Virial 
Coefficient 

The Lennard-Jones potential yields an expansion 
in powers of T~ l/2 [5]. From this key, eq (3) was 
evolved, using notation x=T /T, 

C r =C z 1/2 [l+c£ 3 ][l-exp (1-aT 3 )], 
r =20.615 °K, <7 =1310.5 (cm 3 /g mol) 2 , 
c=2.1486, A=17.4 (cm 3 /g mol) 2 . (3) 

It represents the third virial coefficient within the 
apparent precision of the data. The last factor on 
the right of (3) differs significantly from unity only 
at temperatures below 40 °K. Calculated values of 
C are given in table 3. 
The series 



i=l 



(3a) 



requires at least six terms to give a better representa- 
tion than eq (3). The alternating series obtained by 
least-squares, however, are highly divergent; for 
example, the fifth term for N=6 exceeds the value 
of by factor 240 at 24 °K. 



5. Note on Two-Term Representations 

The behavior of plots of y=B or C as functions of 
x = To/ T suggests an empirical relation of the form 



yly* 



^x fl —x p 



(4) 



with constants m<>. The limitations of (4) having 
been examined in obtaining (2) and (3) above, it is 
reasonable next to examine the relation 



y/yo=exp (to)— exp (kx v ). 



(4a) 



Since the five parameters in (4a) must be found 
tediously by trial, the exploration has been discon- 
tinued with preliminary results in the following 
table. It is concluded that (4a) may be sufficiently 
flexible to offer promise of accurate representations 
of both B and C. 



V 


M 


V 


k 


To 


2/o 


n 


A 


H 


H 


H 


0.6 


110.9 


92.18 


49 


0.20 


C 


v 2 


l 


4 


20 


135 


32 


55 



A=(n— ty-^ly— i/caic 
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6. Appendix 



Table 


I. Virial 


coefficients from extended isotherms of eq (1) 


T, °K 


B, ciu 3 /g mol 


C, cm3/g mol 2 


T, °K 


/>', cmtyg mol 


C, cm3/g mol 2 


24 


-115.3 


5601 


40 


-50.18 


1027 








42 


-46.02 


919 


25 


-106.5 


1952 


44 


-42.54 


1134 


26 


-101.1 


2736 


46 


-39. 25 


1108 


27 


-95.01 


2032 


48 


-36.28 


1103 


28 


-90.36 


2624 








29 


-85. 72 


2561 


50 


-33.49 


1050 








55 


-27.67 


1068 


30 


-80.96 


1725 


60 


-22.93 


1057 


31 


-77.24 


2180 


65 


-18.77 


906 


32 


-73.69 


2642 


70 


-15.30 


809 


33 


-69. 14 


554 








34 


-66. 27 


1549 


75 


-12.57 


854 








80 


-9.96 


767 


35 


-63. 29 


1537 


85 


-7.58 


615 


36 


-60.30 


1402 


90 


-5.57 


559 


37 


-57.51 


1280 


95 


-4.00 


637 


38 


-55.04 


1335 








39 


-52. 66 


1321 


100 


-2.52 


620 



Table II. Polynomial coefficients in Lennard-Jones expan- 
sion for second virial coefficient 



B=J2biT- m - 1)n , (cm'/gmol). (2b) 

,1=1 





3 


4 


5 


6 


61 
&9 
63 

b 5 


4-1. 18579. 102 
-7.32017. 10 2 
-5.21897.103 

0. 2991 


+1.37452.102 
-1.25884.103 
-8. 16662. 102 
-l.13500.10i 

0. 0664 


+1. 45098. 102 
-1.55190. 103 
+3. 04443. 103 
-3. 24136. W 

4-4. 07085. 10* 

0.0615 


+ 1.20054.102 
-3. 70405. 10 2 
-1.78006.10* 
+1.41284. 10 3 
-6.48921.105 
+1.05121. 10« 
0. 0408 
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